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Polarization of the cosmic microwave background, though not yet detected, provides a source 
of information about cosmological parameters complementary to temperature fluctuations. This 
paper provides a complete theoretical treatment of polarization fluctuations. After a discussion of 
the physics of polarization, the Boltzmann equation governing the evolution of the photon density 
matrix is derived from quantum theory and applied to microwave background fluctuations, resulting 
in a complete set of transport equations for the Stokes parameters from both scalar and tensor 
metric perturbations. This approach is equivalent at lowest order in scattering kinematics to classical 
radiative transfer, and provides a general framework for treating the cosmological evolution of density 
matrices. The metric perturbations are treated in the physically appealing longitudinal gauge. 
Expressions for various temperature and polarization correlation functions are derived. Detection 
prospects and theoretical utility of microwave background polarization are briefly discussed. 



I. INTRODUCTION 



Since the initial announcement by the COBE team of the detection of cosmic microwave background temperature 
anisotropies [1], a great deal of experimental activity has resulted in nearly a dozen more anisotropy detections on a 
wide range of angular scales [2,3]. Simultaneously, detailed numerical analysis has sharpened theoretical expectations 
for the temperature anisotropy and its dependence on a variety of cosmological parameters, primarily in the context of 
theories with initial adiabatic perturbations (of which Cold Dark Matter is a special case) [4-7] but also in cosmological 
defect models [8]. While much work remains to be done, the focus of microwave background research has shifted from 
simply detecting anisotropies to creating a detailed picture, both experimental and theoretical, of the anisotropies on 
all angular scale, and using this picture to constrain cosmological models [3]. 

The cosmological information in the microwave background is encoded not only in temperature fluctuations but 
also in its polarization. Since, as discussed below in Sec. II, the source term for generating polarization is fluctuations 
in radiation intensity, generally polarization fluctuations are expected to be somewhat smaller than temperature 
fluctuations. Numerical calculations have confirmed this rough expectation, giving polarization fiuctuations no larger 
than 10% of the temperature anisotropies [9,10]. Greater experimental sensitivity is required to measure polarization 
than temperature fiuctuations, and so far only upper limits have been established. However, polarization fiuctuation 
measurements also have certain experimental advantages over temperature fiuctuation measurements, making first 
detection of polarization within the next few years a reasonable possibility. 

Currently the best polarization limit comes from the Saskatoon experiment [2], with a 95% confidence level upper 
limit of 25 //K between two orthogonal linear polarizations at angular scales of about a degree, corresponding to 
9 X 10~® of the mean temperature. An earlier experiment mapped a large portion of the sky with a 7° beam to 
place upper limits of 6 x 10~^ in linear polarization and 6 x 10"'* in circular polarization for quadrupole and octupole 
variations [11], while a later measurement put limits of around 5 x 10~® on both linear and circular polarization 
at arcminute scales [12]. The COBE satellite made polarized measurements and can in principle achieve a limit of 
around 10~® at angular scales from 7° to quadrupole, although doing so would require reanalysis of the entire data set 
[13]. While the sensitivity of the Saskatoon experiment is very good, it is designed primarily to measure temperature 
fiuctuations. A new experiment planned for 1995 is optimized to measure polarization at 1° and 7° scales and aims 
for a sensitivity at or below 10~® [14]. 

On the theory side, microwave background polarization has been discussed for many years [15-19]. Detailed 
calculations for CDM models have given linear polarization as large as 10% of the temperature anisotropy at medium 
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angular scales, with a strong dependence on the ionization history of the universe [9,10]. Another recent calculation 
has mapped expected correlation patterns between temperature fluctuations and polarization [20]. 

The aim of this paper is a detailed investigation of the theory of microwave background polarization. In contrast to 
previous work employing classical radiative transfer theory [21], the evolution of polarization is derived from a photon 
description; this approach has previously been applied to temperature fluctuations in systematic investigations of 
second-order effects [6,22]. The usual classical Boltzmann equation, adequate for describing temperature fluctuations, 
must be generalized to a density matrix describing the photon polarization state. The formalism for this generalized 
Boltzmann equation was recently developed in the context of neutrino-flavor evolution [23]. Here the appropriate 
equation is derived beginning with the fundamental description of the relevant Compton scattering process; the 
techniques easily generalize to give a Boltzmann equation for any particular density matrix. The advantage of this 
approach is that it treats photons in a general manner, like other particle species described by a Boltzmann equation, 
and can easily be applied to other polarized distributions, e.g. electrons in a magnetic field. It also gives a systematic 
perturbative expansion in the relevant small quantities, and thus provides the framework for an investigation of all 
second-order polarization effects, which may be of particular interest in the case of early reionization for which the 
polarization contribution is the largest [24,25]. This paper also serves as a review which presents many derivations of 
relevant formulas which remain unpublished. 

Section II provides an overview of the physics of polarization and its application to the microwave background. 
A simple calculation demonstrates that only quadrupolar variations in radiation intensity on a scatterer produce 
polarization. The Stokes parameters are defined and their connection to the photon density matrix made explicit. 
Section III derives the general formula for the time evolution of a density matrix in terms of an interaction Hamiltonian. 
This section is rather formal; the relevant result is Eq. (3.12). Section IV specializes this result to the evolution of 
the photon density matrix including Compton scattering. The calculation in this section is straightforward but long; 
the ultimate result is Eq. (4.24), which is equivalent to the usual classical equations of radiative transfer [21]. 

The particular application to a cosmological context begins in Section V, which derives the general relativistic 
Liouville equation for a perturbed Friedmann-Robertson- Walker spacetime. This collisionless part of the Boltzmann 
equation describes the photon geodesies in a homogeneous and isotropic universe with small scalar and tensor metric 
perturbations. The scalar perturbations are described with the physically appealing longitudinal gauge [26] instead 
of the more traditional synchronous gauge. Section VI gives the final evolution equations for the Stokes parameters 
describing the microwave background. The temperature and polarization fiuctuations must be expressed in terms of 
their statistical properties for comparison with experimental results; Section VII derives expressions for the power 
spectra of various correlations and cross-correlations between temperature and polarization. Finally, the concluding 
section briefiy discusses what detailed polarization measurements may eventually reveal about cosmology. 

This paper employs natural units throughout, in which h = c = G = kB = l', Section II uses gaussian units 

for electromagnetic quantities. The metric signature is (-|- — ). Spinor normalizations and other field theory 

conventions in Section IV conform to Mandl and Shaw [27]. 

II. PHYSICS OF POLARIZATION 

This section gives a qualitative overview of the physics of polarization in the context of the microwave background. 
We begin with a review of Stokes parameters, the conventional method for describing polarized light. Then we 
show how polarization can be generated by scattering; application to processes on the last scattering surface predict 
distinctive correlations between hot and cold spots and polarization direction. The equivalent polarization description 
in terms of the photon density matrix is then presented, with the connection to the conventional Stokes parameters 
made explicit. 

A. Review of Stokes Parameters 

Polarized light is conventionally described in terms of the Stokes parameters, which are presented in any op- 
tics text [28]. Consider a nearly monochromatic plane electromagnetic wave propogating in the z-direction; nearly 
monochromatic here means that its frequency components are closely distributed around its mean frequency uiq- The 
components of the wave's electric field vector at a given point in space can be written as 

= a^it) cos [ujot - 6x(t)] , Ey = ay (t) cos [ujot - 6y(t)] . (2.1) 

The requirement that the wave is nearly monochromatic guarantees that the amplitudes and Gy and the phase 
angles and 6y will be slowly varying functions of time relative to the inverse frequency of the wave. If some 
correlation exists between the two components in Eq. (2.1), then the wave is polarized. 
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The Stokes parameters are defined as the following time averages: 



I={ai) + {afj); (2.2a) 

Q={al)-{al); (2.2b) 

U={2a,aycos(0,-0y)); (2.2c) 

V = {2a^ay sm{e^ - Oy)). (2. 2d) 

The parameter I gives the intensity of the radiation which is always positive. The other three parameters define 
the polarization state of the wave and can have either sign. Unpolarized radiation, or "natural light," is described 
hy Q = U = V = 0. One important property of the Stokes parameters is that they are additive for incoherent 
superpositions of waves. The four parameters can be measured with a linear polarizer and a quarter-wave plate; the 
first three can be measured with only a linear polarizer. The V parameter can also be measured as the intensity 
difference between left and right circular polarizations. 

The parameters I and V are physical observables independent of the coordinate system, but Q and U depend on 
the orientation of the x and y axes. If a given wave is described by the parameters Q and U for a certain orientation 
of the coordinate system, then after a rotation of the x — y plane through an angle (f), the same wave is now described 
by the parameters 

Q' = Q cos{2(f)) + Usm{2(f)), 

U'= -Qsm{2(f)) + U cos{2(f)). (2.3) 

From this transformation it is easy to see that the quantity + is invariant under rotation of the axes, and the 
angle 

a = -tan"^ — (2.4) 

2 Q ^ ^ 

transforms to a — cf) under a rotation by (f) and thus defines a constant direction. The physically observable polarization 
vector P is here defined as orthogonal to the direction of wave propogation, having magnitude (Q^ + [/^)^/^ and polar 
angle a. For a wave with linear polarization, the vector P lies along the constant orientation of the electric field. 
Note that since the definition (2.4) is degenerate for a and a + tt, only the orientation of P is defined and not the 
direction. We take the range of a to be — 7r/2 < a < 17/2 with the sign of a the same as the sign of U . While the 
radiation transport equations below are most conveniently formulated in terms of the Stokes parameters, the physical 
interpretation of a polarization pattern is clearest in terms of the observables I, V, and P. 



B. Polarization and the Last Scattering Surface 

In the early universe, at redshifts greater than about z 1100, the baryons, electrons, and photons comprise a 
tightly coupled fiuid. Small metric perturbations induce bulk velocities of the fiuid, and the resulting anisotropies in 
the photon distribution will induce polarization when the photons scatter off charged particles. After recombination, 
the photons freely propagate along geodesies, and any polarization produced before recombination will remain fixed. 
A sufficiently early reionization can of course generate further polarization. 

A simple idealization of the last scattering surface elucidates the process of polarization generation. Consider 
initially unpolarized light which undergoes Thomson scattering at a given point and is then viewed by an observer. If 
the intensity of the light incident on the scattering point is uniform in every direction, then obviously no polarization 
can result; however, if the incident intensity varies with direction then polarization can be generated. Choose the 
z-axis to lie in the direction of the outgoing light, which is described by the Stokes parameters I, Q, U , and V; 
represent the light incident on the scattering point by the intensity I'{9,<j)). Define the polarization vectors for the 
outgoing beam of light so that is perpendicular to the scattering plane and iy is in the scattering plane, and likewise 
with the incoming polarization vectors and e'y (see Fig. 1). Also, instead of dealing with I and Q, it is convenient to 
describe the scattering process in terms of 7^; = (I + Q)/2 and ly = (I — Q)/2. The Thomson scattering cross-section 
for an incident wave with linear polarization e' into a scattered wave with linear polarization e is given by 

where cry is the total Thomson cross section. The incoming wave is unpolarized by assumption, and thus satisfies 
I'^ = I'y = I' /2. The scattered intensities are 
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FIG. 1. Definition of vectors and angles for Tliomson scattering of a ligfit beam or plioton. 



h= ^ K(C • ^yf + Ki^y ■ ^^yf] = ^-^I' cos^ 0. (2.6b) 



Thus the scattered wave has the Stokes parameters 

I=h + Iy = ^I'(l + cos''e), (2.7a) 

Q=L-Iy = ^I'sm^e. (2.7b) 

This calculation gives no information about the U or V parameters. As will be shown later, the V parameter remains 
zero after scattering and will not be considered further [21]. The U parameter can be determined by using Eq. (2.3). 
Simply rotate the outgoing basis vectors in the above calculation by 7r/4 and recalculate Q, which will be equal to U 
in the original coordinate system. The result is = 0. These results can alternatively be obtained from the physical 
description of the polarization state in Rayleigh scattering [21]. Note that Eq. (2.7b) gives the well-known result that 
sunlight from the horizon at midday is linearly polarized parallel to the horizon. 

The total scattered intensities are determined by integrating over all incoming intensities. Note that the outgoing U 
and Q flux from a given incoming direction must always be rotated into a common coordinate system, using Eq. (2.3). 
The result is 

[ dn(i + cos^ e)i'(e,<f>), (2.8a) 

167r J 

Q=^ f dn smHcos{2(f,)I'{9,(f,), (2.8b) 
Idtt J 

U=-^ f dn smHsm{2(f,)I'{9,(f,). (2.8c) 
Idtt J 

The outgoing polarization state depends only on the intensity distribution of the unpolarized incident radiation. 
Expanding the incident intensity in spherical harmonics, 

l'(e,ct,) = J2aimYlm(0,ct,), (2.9) 

Im 



leads to the following expressions for the outgoing Stokes parameters: 



3(Tt 



8 ^ 4 /tt 



-^/^^ aoo + -\l — 020 



(2.10a) 
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FIG. 2. The quadrupolar components of the incident intensity distribution. Any orientation of a quadrupolar distribution 
can be written as the sum of these two distributions. The small arrows indicate the corresponding fluid velocity in a tightly 
coupled fluid. 



Thus scattering generates polarization from initially unpolarized radiation if the radiation intensity at a given point 
as a function of direction has a non-zero component of Y22. 

This particular form for the source of polarization leads to a correlation of the direction of the polarization vector P 
with hot and cold spots on the cosmic microwave sky [20]. Consider a given region on the last scattering surface with 
a spherical mass overdensity; the electron-photon fluid will have a bulk velocity towards the center of the overdense 
region with a velocity gradient away from the center (material further from the center will be falling inwards more 
quickly). In the frame of some particular scattering point away from the center, the fluid velocity in towards the point 
is greater along the radial direction than perpendicular to it, resulting in a quadrupolar radiation intensity variation 
with the largest intensity along the radial direction (see Figure 2). Choose an observation direction at a right angle to 
the radial direction and take this direction to be the polar axis. Then the radiation intensity at the scattering point 
will have a component proportional to cos(2(/) — 2/3) with a positive coefficient, where [3 is the radial direction. The 
scattered Q intensity is proportional to the cos(2(/)) dependence of the incident intensity and the scattered U intensity 
is proportional to the sm{2<j)) piece, by Eqs. (2.9) and (2.10). Thus a, the direction of P in Eq. (2.4), lies along [3, 
the radial direction. For the opposite situation, that of a mass underdensity, all the velocities change sign, so both 
Q and U change sign and the direction of P changes by 7r/2. When the dominant contribution to the temperature 
fluctuations is a gravitational potential difference (Sachs- Wolfe effect [29]), a mass overdensity corresponds to a cold 
spot in the microwave background; in this case cold spots will have radially correlated polarization and hot spots 
tangentially correlated polarization, in agreement with the result of Ref . [20]. For adiabatic acoustic oscillations, the 
density and velocity perturbations are out of phase so no speciflc correlation results. 



The Stokes parameters can be deflned equivalently in terms of a quantum-mechanical description. The polarization 
state space of a photon is spanned by a pair of basis vectors, which we take to be the orthogonal linear polarizations 
|ei) and |e2)- For a photon propogating in the z-direction, the basis states |ei) and |e2) are oriented like the x and y 
axes, respectively. An arbitrary state is given by 




(2.10c) 



C. Photon Description 



(2.11) 
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The quantum-mechanical operators in the linear basis corresponding to each Stokes parameter are given by 



7 = 

Q-- 

U-- 
V-- 



kiX^il- 

kiX^il 

ki>(f2|- 

i\^2){ei\ 



k2>(e2 
k2>(e2 

k2>(ei 

■ «'ki>(e2 



(2.12a) 
(2.12b) 
(2.12c) 
(2.12d) 



The single-particle state expectation values of these operators reproduce the definitions (2.2). For photons in a 
general mixed state defined by a density matrix p, the expectation value for the I Stokes parameter is given by 



(!) =tipi = tr 



Pii 

P'2l 



Pl2 
P'2'2 



Pll + P22 



(2.13) 



and similarly for the other three parameters. These relations thus give the density matrix in the linear polarization 
basis in terms of the Stokes parameters as 



P' 



I + Q 

U + iV 



u ■ 
I- 



iV 

Q 



Q(T3 + Udi + V(T2) 



(2.14) 
(2.15) 



where 1 is the identity matrix and <Ti are the Pauli spin matrices. Thus the density matrix for a system of photons 
contains the same information as the four Stokes parameters, and the time evolution of the density matrix gives the 
time evolution of the system's polarization. 



III. EVOLUTION EQUATION FOR THE NUMBER OPERATOR 

This section considers the quantum number operator for a system of particles and derives its evolution equation, 
including local particle interactions. Taking the expectation value of the operator equation gives the Boltzmann 
equation for the system's density matrix, which is a generalization of the usual classical Boltzmann equation for 
particle occupation numbers (the diagonal elements of the density matrix). The derivation here applies techniques 
previously developed in the context of neutrino mixing [23,30]. 

We adopt second-quantized formalism with creation and annihilation operators for photons and electrons obeying 
the canonical commutation relations 

«!'(p')]= (27r)32/^3(p _ p,)^^^,^ (3^^) 
{K{q), hl,{q')]= (27r)^^^^(q - q')^r.', (3.2) 

where s labels the photon polarization and r labels the electron spin; bold momentum variables represent three- 
momenta while plain momentum variables represent four-momenta. The density operator describing a system of 
photons is given by 

^^^Pii(p')aJ(p')ai(p')> (3-3) 

where pij is the density matrix. The particular operator for which we want the equation of motion is the photon 
number operator 

I?,,(k)=aJ(k)a,(k). (3.4) 
The expectation value of V is proportional to the density matrix, as seen by direct calculation: 

(I?,,(k)) = tr = J ^(p|/5I?,,(k)|p) = (27rf6(0)2k°p,,(k). (3.5) 

The last equality results from repeated application of the commutation relation Eq. (3.2); the infinite delta function 
results from the infinite quantization volume necessary with continuous momentum variables, and cancels out of all 
physical results. 
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The time evolution of the operator T>ij , considered in the Heisenberg picture, is 

^V,,=i[H,V,,] (3.6) 

where H is the full Hamiltonian. We write the Hamiltonian as a sum of the free field piece plus an interaction term: 

H = Ho + Hj (3.7) 

where the interaction piece is a functional of the full fields in the problem. Our goal is to express the right side of 
Eq. (3.6) as a perturbation series in the interaction Hamiltonian Hi. We make the usual assumption of scattering 
theory that in a given interaction the fields begin as free fields and end as other free fields, and the interactions are 
isolated from each other. Consider the evolution of an operator through a single interaction beginning at t = 0: before 
this time, the fields can be taken as free to a good approximation; at t = the interaction Hamiltonian begins to turn 
on, and the interaction finishes at some later time, after which the fields can be taken as free once again. Then the 
time dependence of an arbitrary operator ^ to first order in the interaction Hamiltonian can be expressed as [23] 

^{t)=e{t)+i f dt'[H°{t-t'),e{t)], (3.8) 

Jo 

where C°(t) is the free-field operator with initial condition ^"^(0) = Ci^); ^nd Hj is the interaction Hamiltonian as a 
functional of the free fields. 

Equation (3.8) can be proven as follows. The time derivative of both sides gives 

i^t) = ft^°it) + i mo),m] + i[ di'^^ [HUi - nen)] . (3.9) 

The time derivatives in the first and third term on the right side can be replaced by commutators with H{t) using 
the Heisenberg equation. But these two terms depend only on free fields which are evolved with the free Hamiltonian 
Ho(t) = -ffo(O). Equation (3.9) becomes 

l^it) = i[H°iO),eit)]+i[HoiO),m], (3.10) 

and so to first order in Hi this just gives the Heisenberg equation for the operator ^. 

Now we can express the time evolution oiVij in terms of free field operators. Applying Eq. (3.8) to the commutator 
on the right side of Eq. (3.6) gives 

|l?,,(k) = i [F°(t),I?°.(k)] - j'^dt' [H°i{t-t'), [F°(t),I?°.(k)]] . (3.11) 

The integral on the right side can be cast in a more practical form by making the following physical assumption: the 
duration of each collision (the time interval over which the interaction Hamiltonian is non-negligible, on the order of 
the inverse energy transfer) is small compared to the time scale for variation of the density matrix (on the order of 
the inverse collision frequency). The collision process relevant to the microwave background is Compton scattering off 
electrons, and for the cosmological epoch of interest, the electron density is always low enough for this condition to 
be easily satisfied. Then the time step t in Eq. (3.11) can be chosen large compared with a single collision and small 
compared to the time scale for density matrix evolution. After extending the time integral to infinity and taking the 
expectation value of both sides, we find [23] 

{2^)H{())2k'-p,,{^,) = i{[HUt),v%{^,)])-- J jt{[H°iit),mo),v°^im)- (3-i2) 

Here the integral from zero to infinity has been replaced by an integral over all time; the difference is a principle part 
integral which is second-order in the interaction Hamiltonian and thus ignored. 

Equation (3.12) is the Boltzmann equation for the density matrix pij . The first term on the right side is a forward 
scattering term which is responsible for the MSW effect in a neutrino ensemble [30]; for photons this term is zero, 
as will be shown below. The second term on the right side is the usual collision term. The time integral over the 
free field time dependence enforces energy conservation in each collision. The interaction Hamiltonian will in most 
cases depend on background fields; for example, in the case of the microwave background the Compton scattering 
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collisions are essentially four-point interactions quadratic in both the photon field and the electron field. In principle, 
a second coupled equation for the electron density matrix must be solved simultaneously. However, in many physical 
situations, the background fields may be assumed to have a fixed distribution, generally thermal. In the early universe, 
the electrons maintain a thermal distribution to a very high approximation and the evolution of their density matrix 
becomes trivial. 

The derivation of Eq. (3.12) has been completely general. In the appropriate limit the classical equations of radiative 
transfer are reproduced; the advantage to the current approach is that it provides the same formal framework for 
treating polarized photons as for treating neutrinos or any other particle species governed by the Boltzmann equation. 
It also gives a systematic method for analyzing all higher-order effects. This Boltzmann equation has previously been 
applied to neutrinos interacting through both charged and neutral current processes in supernovae [31]. 



IV. APPLICATION TO COMPTON SCATTERING OF PHOTONS 



In principle, the complete evolution of the cosmic microwave background is determined by Eq. (3.12), generalized 
slightly to include spatial dependence of all quantities. In this work this space dependence will simply be put in 
by hand when taking expectation values and assumed implicitly; more formally it can be included through Wigner 
functions, describing a joint space-momentum distribution [32]. All that remains to be done is substitution of the 
correct interaction Hamiltonian and simplification of the right side. General relativistic terms emerge from the total 
time derivative on the left side; these will be treated in detail in Sec. V. 

Microwave background photons interact with all charged particles. However, the rate of scattering varies with 
the mass of the charged particle as the inverse mass squared; thus it is an excellent approximation to consider only 
Compton scattering off electrons and ignore baryons. This section proceeds with evaluation of the right side of 
Eq. (3.12) for Compton scattering. 



A. Interaction Hamiltonian 



The interaction Hamiltonian density for the fundamental three-point interaction of QED is given by [27] 

Wqed(«) = -e : il>(x)^(x)il>(x) : (4.1) 

where ip is the electron field operator, A'' is the photon field operator, a slash indicates contraction with 7^, and the 
colons signify normal ordering of the operator product. The interaction Hamiltonian is the density integrated over all 
of space: 



HqEB(t) = J (fimqMx). (4.2) 
The scattering matrix describing all scattering processes in QED is given in terms of the interaction Hamiltonian by 

S = J2S("'> = J2--^ d\i...d\„T{nqEB{xi)...nqEB{xn)} (4.3) 

n n * '-' 



n=0 n=0 



where T signifies a time-ordered product. The nth term in the series represents all scattering processes with n 
interaction vertices. Compton scattering is thus contained in the n = 2 term of the scattering matrix. Comparing the 
n = 2 term with the n = I term gives the interaction Hamiltonian for second-order scattering processes: 

5'(')=--/ dt dt'T{H(^EB{t)H(^EB{t')} 

= -i / dtH(^\t). (4.4) 



Using Wick's theorem to simplify the time-ordered product and denoting the piece of iJ^^) describing Compton 
scattering by Hi yields 

Hi(t) = e^ I rft' /rfVv;-(a;)7«5'F(a;-a;')7'' U^(*)-4j(«') + -4«(«')-4j(«)l t/'+(«') (4.5) 
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where Sp is the Feynman propagator for the electron, and ip'^ (V"") ^nd A'^ are linear in absorption (creation) 

operators of electrons and photons respectively. Fourier transforms of the fields and propagator are defined using the 
following conventions: 



^"("^^ = ' (^^^E + , (4.6a) 



(27r)3 k° 

Sf(x) = I ,fl,, ^^".^ e-^'\ (4.6d) 
J (/tt)^ — + 

where Ur(k) is a spinor solution to the Dirac equation with spin index r = 1,2 and e'^ik) are photon polarization 
four-vectors, chosen to be real, with index s = 1,2 labeling the physical transverse polarizations of the photon. The 
summation convention over repeated spin and polarization indices is always implied. The Fourier-space interaction 
Hamiltonian is obtained by substituting Eqs. (4.6) into Eq. (4.5). The distributional identity 

i'^xe'^'^ = (2TTy6\k) (4.7) 

allows trivial integration over the four-momentum of the electron propagator. The resulting interaction Hamiltonian 



IS 



H°j(t)= I rfqrfq'rfprfp'(27r)3,53(q'-Fp'-q-p)exp [it(q'° + p'° - q° - p°)\ 

X \bliq')al{p') {Ml + M2) a.ipMq)] , (4.8) 



M=Mi+M2, (4.9a) 

M,(q'r',p'.s',qr,p.s) ^ .^ -Aq^ A,p')(J> + i + rn)^ MM^) ^ (4 g^) 

2(p ■ q) 



// / / x_ 2«r'(9')^.'(p')(^-/ + "^)^.(p)«r(g) 



M2{q'r',p's',qr,ps) = -e' - v. v. /v. . /..v^z-v./ ^ ^^ g^^ 

2(p' ■ q) 



with the abbreviations 



for electrons and photons respectively. All of the operators in Eq. (4.8) are free-field operators, so this is the proper 
expression to substitute into the left side of Eq. (3.12). 
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B. Forward Scattering Term 



We now proceed to evaluate the first term on the left side of Eq. (3.12). First we display operator expectation values 
needed here and in the following subsection, using operator definitions and the commutation relations, Eq. (3.2): 



(aia2 • • • &1&2 • • •) = («i«2 • • •) (&1&2 • • •) 
(a|„(p')«n(p)) = 2p°(2iTf6^(p - p')p™n(p) 



(4.11a) 
(4.11b) 



m 2 



(4.11c) 



(&):;(9'l)&r,(?l)&:-(?2)&r,(?2))= ^(27r)«^^(qi - q'i)^^(q2 - q'2)^.,./^.,.,' f na(qi)na(q2) 



I1I2 



2 1 2 



1 - 2"«(qi) 



{al,y,)a,,ip^)al,ip',)a,,ip2))= ipy,i27rf6'ipi - p'i)^^(p2 - p'2)p.,.;(pi)p.,.^(P2) 



(4. lid) 



(4. lie) 



The last relationship neglects the correlation term between all four operators when p'l = pi = p'^ = p'l- The 
expectation values for electron operators assumes a particular form for the electron density matrix appropriate to 
thermal equilibrium, with equal populations in each spin state and no correlations between the states; ne{q) represents 
the number density of electrons of momentum q per unit volume. This assumed form for the electron density matrix 
will not be valid if substantial magnetic fields are present. 

Using the definitions (3.4) and (4.8) and the commutation relations (3.2) the commutator in the forward scattering 
term becomes 

[F^°(0),I?°.(k)] = y" rfqrfq'rfprfp'(27r)3,53(q' + p'-q-p)(A^l+A^2) 

X \hl\<i)K{<l)a\\f)a^{\z)2p\2i,fh,b''{j> - k) 

-hl\<i)K{<i)a\{\z)a,{pY2p'\2i,fb^A''{j>' - k)] . (4.12) 
On using the above expectation values, it follows that 

i([F°(0),I?°.(k)]): 



T 



c^q ((^»s/Os'j(k) - (5js'/9is(k)) 



h ■ q 



(4.13) 



where the integrals have been performed with the delta functions. All of the terms involving 1/: cancel out on using 
the gamma-matrix identity jfilj^ = 2 A ■ B — Ij^jf^ and the polarization vector properties k ■ ei{k) = and Ei ■ Ej = —bij . 
For the remaining terms we use the identity 

^r{q)^,,{4 + m)^^,Ur{q)= Ur{q){2q ■ Ss - + m^J^^,Ur{q) 

= 2q ■ SsUr{q)^s''^r{q) 



— {q ■Ss)(q -Ss') 
m 



(4.14) 



where the second equality follows from the Dirac equation and the third equality uses the Gordon identity. Thus we 
have 



i([F°(0),I?°.(k)]) = 

and the forward scattering term does not contribute to the photon density matrix evolution. 



(4.15) 
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C. Scattering Term 



The scattering term is considerably more cumbersome to evaluate, being quadratic in the interaction Hamiltonian. 
After substituting the expressions for Hi and Vij and taking the expectation value, the scattering term reads 

i dt{[H°{t),[H°{0),V°^{k)]]) = j rfqirfq;rfpirfp;rfq2rfq'2rfp2rfp'2 

X (27r)7^3(q, ^ p, _ _ p^)^4(^, ^ _ _ p2)A^(l)A^(2) 

-rf(27r)3^,,/3(pi - k) (gy^(g2)) («!;(K)%(fc)«!^(P2)«.,(P2)) 

-p'i°(27rf ^,.'/3(p'i - k) (gy^(g2)&J;(9'i)^(9i)) («!^(P2)«.,(P2)«IW«.,(Pi)) 
Vi°(27rf^,,;^^(p;-k)(&t,(5;)6^^(5^)6t,(5^)6^^(5^)^^„t(^)^ (4,16) 

The energy delta function comes from the time integral on using Eq. (4.7). The arguments of the matrix element 
indicates the subscript to be attached to all dependent variables in Eq. (4.9), and of course summation over all spin 
and polarization indices is implied. 

Substitution of the expectation values Eq. (4.11) into the above expression and performing the integrals over q2, 
P2, q2 and p2 using the various delta functions yields 

- / rft([F°(t),[F°(0),I?°.(k)]]) 

^ J -co 

= i(27r)3,53(0) j d(id(i'dp'{27rf6^{q'+p'-q-k)M{q'r',p's[,qr,ksi)M{qr,ks'2,q'r',p's2) 

+ ^{2Trf6^{0) j d(id(i'dp{2Trf6^{q' + k-q- p)M{q'r', ks[, qr, psi)M{qr, ps'^, q'r', ks2) 

X [ne(q')'^is'/sjs^fts,(k) - ne{q)6js[6is.^Psl^sSv)\- (4-17) 

The subscript "1" on all momentum variables has been dropped for notational simplicity. All terms quadratic in the 
electron phase-space density have been dropped since for all cosmological scenarios this number is negligible compared 
to unity; all terms quadratic in the photon density matrix cancel exactly. By relabeling the integration variables and 
spin indices (implicitly summed over) in the second integral, Eq. (4.17) reduces to 

- / rft([F°(t),[F°(0),I?°.(k)]]) 

= i(27r)3,53(0) j d(id(i'dp'{27rf6^{q'+p'-q-k)M{q'r',p's[,qr,ksi)M{qr,ks'^,q'r',p's2) 

X [ne{(i)6,^,>^ {6i,^p,>^j{k) + 6j,>^pi,^{k)) - 2ne(q')'^is/is^Ps;s,(p')] • (4-18) 

This equation is an essentially exact expression for the collision term in the case of the microwave background: the 
approximation that the duration of the Compton scattering be small compared to the time between scatterings is 
eminently satisfied for any cosmological scenario, and assuming the electrons to be unpolarized is essentially exact 
unless magnetic fields become important at some epoch. 

Evaluating the matrix elements and performing the integrals in Eq. (4.18) is a straightforward process. This paper 
is concerned with the first-order perturbations away from a perfectly homogeneous and isotropic universe, and the 
scattering term will be explicitly calculated to first order. Evaluating the matrix elements involves standard techniques 
of quantum field theory and yields the familiar Compton cross-section to lowest order: 
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^M{q'r' ,ps[,qr, ksi)M{qr, ksij, qr',ps'^) 



2e^ 



q ■ p q ■ 1^ / 



+2(£,^(k) • £s;(p)£s,(p) •£s^(k) +£,^(k) • £s,(p)£s;(p) •£s^(k) - (5,^,^(5,'^, J 



(4.19) 



The following subsection then obtains the general Boltzmann equation for the photon density matrix to first order in 
terms of the photon energy and polarization vectors. 

D. Scattering Term to First Order 

Now we proceed to evaluate Eq. (4.18) to lowest order in scattering kinematics. After substituting the matrix 
element Eq. (4.19) into Eq. (4.18), the Boltzmann equation (3.12), now explicitly including spatial dependence, 
becomes 



,,P8j(x,k)— 
at Ibm^k 



dq dp- 



-{27r)6{E{(i + k - p) + p - Eiq) - k) 



E(q+k-p) 

X (ne(x,q)6,^,^^ {6i,^p,y(x,k) + (5^,^^,^ (x, k)) - 2ne(x,q')6i,^6j,^^p,^^,Jx,p') 
q ■ p q ■ k 



, + [ss[ip) •£s^(k)£s^(k) • £,^(p) - £,^(k) •£s,(p)£s;(p) • £s^(k) + 

q ■ K q ■ p J ^ 

+2(£.'(p) (k)£,- (k) •£ (p) + £ (k) (p)£,'(p) •£.'(k) - 8,1 8_ 



1^2 12 



(4.20) 



where -E'(q) = (q^ + rn^y/^ is the energy of an electron with momentum q. The electrons are described by an 
unpolarized thermal Maxwell-Boltzmann distribution: 



ne(x, q) = ne(x) 



\mTe 



3/2 



exp 



(q — mv(x))^ 



2mTe 



(4.21) 



with Tg the electron temperature and v(x) the electron bulk velocity. Useful integrals of the electron distribution are 



;"e(x,q) = ne(x). 



(4.22a) 



-qine(x,q) = mvi(x)ne(x). 



(4.22b) 



For relevant cosmological situations, the kinetic energies of the electrons and photons are negligible compared to the 
electron mass, implying that the energy transfer in a Compton scattering event is small compared to the characteristic 
photon energy: p <^ m, q <^ m (with obvious abbreviations p = p° = |p| and q = |q|). Furthermore, if the electron 
and photon temperatures are comparable, p <^ q. We expand the various functions in Eq. (4.20) in terms oi p/q and 
q/m, using the following asymptotic expansions: 



^(q + Q) 



q2 , q Q , Q2 



2m? 



(4.23a) 



"e(q + Q) ~ "e(q) 



Q-(q-mv) 



2mTe 



(4.23b) 
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8{k-p + ^(q) - ^(q + k - p)) ~ 8{k -p)+ (k v)-^db{k p) ^ ^ ^ ^ ^ 

m op 

where in the last expression the derivative of the delta functional is defined through integration by parts [33]. Writing 
out the polarization sums explicitly yields the equation 



e4n,(x) , f dn 



6(k - p) + (k - p) •v(x) 



d6{k — p) 



dp 



x|-2 + ^ ) /9ij(x,k) + 4p • £i(k)p •£i(k)/9ij(x,k) + 4p •£i(k)p • £2(k) /92i (x, k) 
^ + ^ - 2 ) % (pii(x,p) + P22(x,p) 

f + ~ ) •£i(p)£i(k) •£2(p)-£i(k) •£2(p)£i(k) • £i(p) ) ( pi2(x, p) - p2i(x, p) 



+2(^£i(k) •£i(p)£j(k) •£2(p) + £i(k) •£2(p)£i(k) •£i(p)j (^pi2(x,p) +p2i(x,p)^ 

+4£i(k) •£i(p)£j(k) •£i(p)pii(x,p)+4£i(k) •£2(p)£i(k) • £2 (p) /922 (x, p) | (4.24) 

Here the photon momentum integral has been rewritten as an energy integral and an angular integral over the 
momentum direction. This is the basic equation describing the evolution of the photon density matrix to first order 
in the kinematic variables. By rewriting the momentum and polarization vectors in a spherical coordinate basis and 
incorporating the velocity-dependent term into the left-hand side, the equation becomes equivalent to Chandrasekhar's 
radiative transfer formalism; cf. Chapter 1, Eq. (212) of Ref . [21]. Before performing the final angular integrals, we 
must consider the left side of the equation for the particular space-time geometry in which we are interested, which 
determines the azimuthal dependence of p. The left side of the equation will be analyzed in the following section, and 
then in Section VI we perform the remaining momentum integrals to complete the evaluation of the right side. 

V. THE GENERAL-RELATIVISTIC LIOUVILLE EQUATION 

The last section has analyzed the right side of the Boltzmann equation, Eq. (3.12); we now turn to the left side, 
describing the propogation of photons in the background space-time. The Boltzmann equation with no collision term 
on the right side is the Liouville equation, describing the evolution of a collisionless system's phase space distribution. 
Writing the equation has already assumed definition of a set of space-like hypersurfaces; that is, the equation contains 
an explicit time derivative. The background space-time here will be the canonical Friedmann-Robertson- Walker 
(FRW) space-time. In this paper only the fiat case will be considered; techniques pertaining to open universes have 
also been extensively developed [34]. Scalar and tensor metric perturbations are added to the fiat background space- 
time; we neglect vector perturbations, which kinematically decay and are unimportant unless a continual source of 
vector perturbations, such as topological defects, exists. The metric we consider is 

(/oo = l + 2$(x,t), ffoi = 0, (5.1) 



Qij = -a'{t)[{l - 2*(x,t))% +%(x,t)]. (5.2) 

In this section, Greek subscripts refer to space-time indices running from to 3 while Roman subscripts refer to 
spatial indices running from 1 to 3. The function a{t) is the usual cosmological scale factor. The scalar perturbations, 
defined in the longitudinal gauge, are given by the two scalar function $ and ^ [26]. The metric perturbations hij 
are defined in the transverse-traceless gauge and are subject to the constraints 

hi' = Q, d^hij=Q. (5.3) 

With the perturbations defined in this way, no residual gauge freedom remains, in contrast to the more conventional 
synchronous gauge condition. A second advantage of the present definition is that in the Newtonian limit, the metric 
perturbation $ simply corresponds to the Newtonian potential. The inverse metric to first order in the perturbations 
is 
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g°° = l-2'^{i,,t), g°' = 0, (5.4) 

g'' = -^[(l + 2*(x,t))^'^- -/j'^-(x,t)]. (5.5) 

In this section we will consider the Liouville equation to first order in the metric perturbations; second-order treatments 
have been undertaken in Ref. [35]. 

Photons are described by space-time coordinate and four-momentum . Our coordinate system has x° = t; 
the photon momentum satisfies 



The photons obey the geodesic equation 

£x^' „ dx" dxP , , 

where A is an affine parameter along the photon geodesic, which may be defined so that dtjd\ = dx°/d\ = thus 
dx'/dX = k' using Eq. (5.6). Therefore, using the definition of the Christoffel symbol F, 

dk'^^..[ldgai_dg^\k^ 

dt ^ \2 dx" dxp ) k° ^ ' 

with the geodesic condition k^kfj_ = Q . 

Liouville's equation for any phase-space distribution function / is 

di^d^dx^ d^dk^^di d^k^ d^dk^^ 

dt dxi' dt dki' dt dt dxi k° dki' dt ' ' ' 

We change to a convenient choice of momentum variables, the photon energy in a local orthonormal frame k = 
(—k'kiY/^ and the unit vector k\ A local observer at rest with respect to the cosmic coordinate system will measure 
the photon to have energy k; this is the energy appearing in the collision term on the right side of the Boltzmann 
Equation. The wave vector is thus given by 

F = ifcF (^1 + *- , (5.11) 

and the Liouville equation in the new variables is 

^^di dik^ didk didk^^ 

dt dt dxi k° dk dt dki dt ' ' ' 

Now it is a straightforward matter to calculate the derivatives dk / dt and dk'/dt. Substituting k° = k(l — $) and 
Eq. (5.11) into the zero component of Eq. (5.9) gives to first order 



dk_ 
'dt 



d 5* . k' 5$ . 1 f , dhi 



a dt a dx^ 2 dt 



(5.13) 



The derivative of A;' is most easily computed by differentiating Eq. (5.11) and equating this result with Eq. (5.9); it is 
straightforward to show that dk' / dt has no lowest-order terms. Physically, this is because geodesies are straight lines 
in the unperturbed metric. Since df /dk' is itself linear in the metric perturbations, the final term of the Liouville 
equation drops out to first order. Expanding the distribution function as 

/(x, k, k, t) = f(°\k, t) + /(i)(x, k, k, t) (5.14) 
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leads to the zeroth-order equation 



5/(0) a,5/W ^ 



dt 



a dk 



(5.15) 



whose solution is simply f^'^\k, t) = f^'^\ka), which is just the uniform redshift of the spectrum with cosmic expansion. 
The first-order Liouville equation is 



5/(1) 5/(1) y 



dt 



dx^ a a dk 



a^dj^ d_f^^ 



dk 



dt a dx' 2 dt 



0. 



(5.16) 



Note that k in this equation is the physical, not the comoving, photon wave number. 

The terms in the Liouville equation depending directly on the metric perturbations determine the form of the 
directional dependence of the distribution function. A Fourier transform over the spatial dependence of the equation 
gives for the Boltzmann equation 



^/(i)(K, fc, k)+^(K . k)/(i)(K, fc, ic) - ^fc A/(i)(K, fc, ic) 



I 

-I 

a 

df(°Hk) 



dk 



^<I>(K)-^*(K) + l(K.k)<I>(K) + iF 



dt 



%(K) 



C(K,k,k) 



(5.17) 



where C represents the collision term on the right side. For scalar perturbations, the previous section shows the 
right side contains source terms proportional to k • v, where v is the local velocity of the electrons. But for scalar 
perturbations, v cx K [26], so if we choose spherical coordinates for k with axis in the K direction, then /("'^) is 
manifestly independent of the azimuthal angle (f); in other words. 



/(i)(K,fc,k) = /(i)(K,fc, 



(5.18) 



for scalar perturbations. 

Tensor perturbations do depend on (f), but in a simple manner. We neglect any electron velocity v arising from 
tensor perturbations as corrections to the scalar-induced velocity. The (f) dependence of the distribution function is 
determined by the perturbation term, which can be written as 



FF^%(K,t) 



h+(K,t)e±(K) + h-(K,t)efAK) 



(5.19) 



where efj and efj are polarization tensors for the plus and cross polarizations of the gravity wave. Again, choose 
spherical coordinates with the z-axis pointing in the direction of K. In this coordinate system the polarization tensors 



are given by ej^ 



yy 



1 and e 



with the polarization tensors gives 



xy 



yx 



1, with the other components zero. Contraction of the unit vectors 



*2 2 *2*2 *2 

sin cos (j> — sin sin (j) = sin 0cos2(j>, 



k k-' efj = 2 sin^ cos (f> sin (f> = sin^ sin 2(f>. 
Therefore, for a given plane wave component of a metric tensor perturbation, 

/(i)(K,fc,k) = /(i)(K,fc,(?) cos2^ 
for the plus polarization of the gravity wave and 

/(i)(K,fc,k) = /(i)(K,fc,(?)sm2,/) 

for the cross polarization. 



(5.20) 
(5.21a) 
(5.21b) 



VI. COMPLETE POLARIZATION EQUATIONS 

Now we have assembled all the ingredients for deriving the final polarization evolution equations: Eq. (4.24) with 
the perturbation expansion for pij given by Eq. (5.14) and the angular dependence of pij given by Eqs. (5.18) and 
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(5.21). Note the unperturbed photon density matrix satisfies p'f^^ = and p'f2 = P21 — since it represents a 
uniform unpolarized blackbody spectrum. The lowest-order term on the right side of Eq. (4.24) is zero; the resulting 
equation for /O^*^) gives the uniform shift of the photon spectrum with scale factor in an expanding universe, Eq. (5.15). 
The first order term gives 



4 - 
e We 



dpp 



6(k — p) 



-2i£+i)4'(K.., 



+4p • £,(k) p • £i(k) pg)(K, k) + 4p • £,(k) p • £2(k) p^ffiK, k) 

^ + [e,(k) ■ £l(p)£,(k) . £2(P) - £,(k) • £2(p)£,(k) • £i(p)) (/iV(K,p) - p^V(K,p) 



+2 ( (k) • £1 (p) e, (k) • £2 (P) + (k) • £2 (P) e, (k) • £1 (p) ) ( (K , p) + pI^V (K , p) 



+4£,(k) • £i(p) £, (k) • £i(p) /iV(K, p) + 4£,(k) • £2(p) £j(k) • £2(p) p^V(K, P) 



-(k-p).v(K) 



d6{k — p) 



dp 



+4p • £,(k) p • £i(k) flf^ik) + 4p • £,(k) p • £2(k) pg)(fc) 

+4£,(k) • £i(p) £,(k) • £i(p) pf^ip) + 4£,(k) • £2(p) £, (k) • £2(p) fl^ {p) 



(6.1) 



As in the previous section, K is the Fourier conjugate of x, and is the mean electron density which is constant to 
lowest order. The remainder of this section evaluates the remaining angular integrals in this expression and converts 
the equations for the density matrix elements to equations for the brightness of each Stokes parameter. 

To evaluate the angular integrals most conveniently, choose the z-axis of the spherical coordinate system to coincide 
with K, independently for each Fourier mode. Note that on transforming back to real-space coordinates, care must 
be taken because the density matrix is not invariant under a change of basis (see Sec. II). The basis for the photon 
direction and polarization vectors is taken to be 



sin 9 cos 4> £ij;(k) = cos 9 cos 4> £2j;(k) = — sin ( 
sin sine/) £ij,(k) = cos sin (/) £2j, (k) = cos (/) 



cos ( 



£i.(k) 



■ sm ( 



£2.(k) = 0. 



The same definition is used for p and its associated polarization vectors, with 9^9' and (f) 
integral in Eq. (6.1) is over 9' and (f)' , and the various dot products are given by 



(6.2) 

6'. The angular 



p' 


•£l(k) = 


sin 9' cos 9 cos((/)' — 


<j)) — cos 9' sin 


p' 


•£2(k) = 


sin 9' sin((/)' — (f)), 




£l(k) 


£l(p) = 


cos 9 cos 9' cos((/)' — 


<j)) -\- sin 9 sin t 


£l(k) 


£2(p) = 


— cos 9 sin((/)' — (/)), 




£2(k) 


£l(p) = 


cos 9' sin((/)' — (/)), 




£2(k) 


£2(p) = 


cos((/)' — (f)). 





Two additional convenient abbreviations are = v • p = cos 9' and = v • k 
straightforward, resulting in expressions like 



(6.3) 

Now the angular integrals are 



4 - 



at 167rm^ 



2p' 



-2p' 



(0) 



dk ' ""^"'^ dk 



2kvp^ 



(6.4) 



for scalar perturbations, where dependence on the Fourier mode K is implicit. In solving the evolution equations, it 
is convenient to split the density-matrix perturbation into two parts, one due to the scalar metric perturbations and 
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one due to the tensor metric perturbation. In making this split, the bulk velocity v is entirely attributed to the scalar 
perturbations. 

For the final set of evolution equations, we change variables to comoving wave number q = ka, and convert the 
density matrix elements to Stokes parameter brightness perturbations: 



A}(K,g,//)E 
A'Q(K,g,//)E 
A^(K,g,//)E 
AV(K,g,//)E 



9 dp[°^(q) 



4 dq 
Q dp[°^(q) 



Pn(^,(l,p) + P22(^y(lyP) 



4 dq 
1 (9) 



Pn''(K,9,/^) -P22''(K,9,/^) 



4 dq 

9 dp^°^{q) 
4 dq 



Pi2''(K,9,/^) + P2i''(K,9,/^) 



Pi2''(K,9,/^) - P2i''(K,9,/^) 



(6.5a) 
(6.5b) 
(6.5c) 
(6.5d) 



where the superscript i stands for s,+, or x, representing the three types of possible metric perturbations: scalar and 
two polarizations of tensor. For linear perturbations considered here, the distribution function undergoes no spectral 
distortions and the perturbations are blackbody; in this case, A//4 is just the temperature fiuctuation 6T/T0 with 
To the mean temperature. We also define moments of these variables: 



1 dfi' 



Pj{l,')A\{q,l,') 



(6.6) 



where Pi is the Legendre polynomial of order /. Note these moments are sometimes defined differently [36,37] as 

Mp) = E(-0'A,T,. 



For scalar perturbations, the brightness is governed by the set of equations [17] 



dt a 



— iA//$ 

a 



dt 
dA 



A} - Alo + iv^ - -P2{p){A], + A^2 - A^o) 



Q 



dt 
dAi 

dt 
dAi 

dt 



-iK^Ag-- 
iK/iAlj-- 



-iKiiAy-- 
a 



-(TTfleAlj, 



A^ + 77(1 - ^2(//))(Al2 + A^2 - A^o) 



A^ - -pAl., 



(6.7a) 
(6.7b) 
(6.7c) 
(6.7d) 



The evolution of the brightness thus does not depend on the direction of K, only on its magnitude; the brightness 
depends on the direction of K only through initial conditions, which factor out of the linear evolution equations. The 
equations for U and V have no source terms, so for each K the evolution leaves U = and V = 0. The coordinate 
dependence of Eq. (2.3) gives a non-zero U on transforming back to x space, but V remains zero. 

For tensor perturbations, the evolution equations take their simplest form after the coordinate transformation [38] 



A+= (1 - //2)cos(2^A)A+, Af = (I- iJ^)sm(2<j,)Af, (6.8a) 

A+= (l + //2)cos(2,/))A+, Ag = (l + //2)sin(2,/))A5, (6.8b) 

A+= -2//sin(2<?i)A+, A^ = 2// cos(2<?i)A^. (6.8c) 

The (f) dependence is determined by Eqs. (5.21a) and (5.21b), and the /j, dependence is chosen to simplify the final 
equations. After this change of variables, the brightness equations become [40] 

/)A+ 1 ~ dh+ f~ ~ \ 

^ + iz-A>A+ -2^ = -cyTUe (a+ + A+) , (6.9a) 
5A+ 1 

+ -'^^>^Q= --^Tne (a+ - A+) , (6.9b) 
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A+ = A+, (6.9c) 
+ 



dA 



- ^ + -iKfiA+= -aTneA+, (6.9d) 
ot a 



A+ = -— A+ + -A+ A+ + — A+ + -A+ + -A+ (6.10) 



The definition of /}+ is given in Eq. (5.19). The x tensor perturbation gives same equations. Again, V = since it 
has no source term. 

For a given cosmological scenario, which determines the metric perturbations, Eqs. (6.7) and (6.9) must be evolved 
numerically. This can be done efficiently by expanding the /j, dependence of the brightness functions in terms of 
Legendre polynomials, Eq. (6.6), giving a large set of coupled ordinary differential equations [17]. Several detailed 
codes to calculate temperature fluctuations have been implemented using this scheme [4-7]. 

VII. POWER SPECTRA 

Numerical solution of the above transport equations gives the Fourier space brightness functions A|;(A') and Aqi(K) 

for scalar perturbations and A|;(A') and Afji(K) for tensor perturbations, where e represents the two gravity wave 
polarization states + and x . The temperature fluctuations in real space are then 



l + lEE(2'+l)^'(cos^?')e 



To 4 
" K ; 

X [a|;(K) + sin^ 0' cos 2^A'A+ (K) + sin^ 0' sin 2<j,'Afj(K)^ , (7.1) 

where To is the mean temperature of the microwave background, and {0' ,4>') represents the same direction as {0,<j)) 
except in the coordinate system deflned by the K direction. 

The polarization is more complicated, because for each K mode the coordinate system in the direction {0, <j)) has 
a different orientation; when the Q and U brightnesses are summed up, the axes must be rotated to the orientation 
in the x coordinate system, using Eq. (2.3). To determine this rotation angle for each K mode, refer to Fig. 3; the 
needed angle is labelled the angle between the vectors and 0' . Let the direction of K be denoted by (0k,4'k)- 
On the unit sphere, the lengths of the sides of the spherical triangle ABC are just the angles they subtend. The angle 
0' is given by the law of cosines, 

cos 0' = cos cos 0K + sin sin 0k cos((/)if — (f)), (7.2) 

and the rotation angle is 

sin^' = sin (^if sin((/)if — </') CSC 6*'. (7.3) 
Then the Q and U brightnesses are given by 

^%^=iEE(2'+l)^Kcos^^')e'-- 
° K ; 

X [(^A^;(K) + (1 + cos2 6'')cos2<?i'A+;(K) + (1 + cos^ 6*') sin 2<?i'Ag;(K)) cos 2^' 

+ (-2 cos 0' sin 2(f)'A+i{K) + 2 cos 0' cos 2<?i'Ag;(K)) sin 2^'] , (7.4) 

and U the same except for the replacements cos^' — sin^' and sin^' cos^'. From these two quantities, the 
polarization vector P follows from Eq. (2.4) as 



P(x, (?,,/)) 1 



Ar 



Ap(Ap + Aq) + (t,-——JAp(Ap - Aq) 



To "4^2 L ' ' -|A 



(7.5) 



Ap = JAl + Al. (7.6) 
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FIG. 3. Angles and directions for determining the orientation of different spherical coordinate bases at a given point. 



The predictions of a given cosmological scenario are only statistical. The traditional statistical measure of temper- 
ature fluctuations is the angular power spectrum C'(6), defined by 

+ qi-q2 = COS^?, (7.7) 

where the angle brackets represent an ensemble average over initial conditions; this average can be replaced in calcu- 
lations with an average over space, assuming ergodicity. Many analogous correlation functions for various polarization 
variables are possible: 

(7.8a) 
(7.8b) 
(7.8c) 
(7.8d) 

roughly in increasing order of difficulty to measure. The various combinations of Stokes parameters can also be used 
to form correlations: 

C^^(^)^(^^), (7.9a) 
t^^-(^)-(^^), (7.9b) 

and so on. These correlation functions depend on the orientation of the axes used to define the Stokes parameters. 
Their advantage is that they are easier to calculate than Eqs. (7.8) and are easier to measure when signal to noise 
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is low; their disadvantage is that their physical interpretation is less easily visualized. Correlation functions are 
commonly characterized by the coefficients C; of an expansion in Legendre polynomials: 



c^^{e) = Y.'-^crPi{'^o.e) (7.10) 

— ' 47r 

1=2 

and likewise for the others. The / = 1 term, indistinguishable from the Doppler shift from proper motion with respect 
to the rest frame of the microwave background, is ignored. Note this conventional normalization of C; unfortunately 
differs by a factor of At: from the conventional normalization of the brightness moments, Eq. (6.6). 



A. Temperature Correlation Functions 

The temperature correlation function C'^'^{6) can be evaluated exactly in terms of the brightness moments. First, 
since the equations for A(K) depend only on |K|, separate out the angular dependence by writing the fluctuations as 

A},(K)=A},(A%-(K)e'*-W, (7.11) 

where i stands for s, +, or x. The functions Q and 6i are random variables set by the initial conditions: Q is real and 
positive, normalized such that J dilxCii^) = 47r, and the phase 6i is a real number between and 27r. For gaussian 
initial conditions, Q is drawn from a normal distribution and 6i is uniformly distributed, independently for each K 
mode and each type of perturbation. 

Substituting Eq. (7.1) into Eq. (7.7) provides the starting point for evaluating the temperature correlation functions. 
The average value can be replaced by J d^x, where V is the sample volume, and in the limit of a large volume 
the sums over K vectors can be replaced by integrals: V/i^ir)^ f rf^K. The x integral gives a delta functional 

which eliminates one of the K integrals. For scalar perturbations, the result is 



mi)T(q2)\ V 



- I rf3K^(2/i + l)(2/2 + l)P,Acose[)P,,(cose',)A}l(K)A},jK). (7.12) 



To To I 128.. 

Then expressing each Legendre polynomial in terms of spherical harmonics using the spherical harmonic addition 
formula and performing the angular piece of the K integral using the decomposition in Eq. (7.11) gives the familiar 
formula 



rprp _ V 



- / A-2rfA-|Al,(A')p. (7.13) 



For tensor perturbations, the added angular dependence makes the calculation somewhat more involved; in this 
case, Eq. (7.12) becomes 



r(qi)m2)\_ y 



To To I 1287r3 ., 

'1 '2 



y" rf3K^(2/i + 1)(2/2 + i)i^;i(cos e'^)Pi,{zo^e'^)^Y^\e'^)^Y^\e'^) 



X [cos 2<i>\ cos 2,/-'2 A+ * ( A') A+ ^ ( A)C+ (K)^ + sm 2<i>\ sm 2,/-'2 A]<,* ( A) A]<,^ ( A)Cx (K)^] . (7.14) 



The cross-terms between the two types of tensor perturbations (and between the tensor and scalar perturbations) 
cancel because the integral over K contains random phase factors, e.g. exp(i(5_|_ (K) — ib^, (K)). Rearranging the terms 
in brackets yields 



/ r(qi)nq2) \_ y 

\ To To I 128.. 

1 1 I 



xi [cos(2,/-; + 2^',) (a+:(A)A+JA)C+(K)2 - A]<,:(A)A]< jA)Cx (K)^) 

+ cos(2,/-; - 2^',) (a+:(A)A+JA)C+(K)2 + A]<,:(A)A]<jA)Cx(K)2)] . (7.15) 

Further progress can be made with the mild and physically reasonable assumption that Aj;(A') = Aj;(A'), in other 
words that the power spectra for the two polarizations of gravity waves are the same. Then the two terms on the 
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second line of the above equations give equal integrals and cancel. The remaining trigonometric functions can be 
written in a form independent of 4>' [41]: 

T{(li)T{(i2) 

To To 



-|^^(2/i + l)(2/2 + 1) / dKK'Aj*{K)AjjK) 



X / dnMl,[)PM) [2(qi • q2 - - (1 - - 1^?)] 



(7.16) 



where as in the scalar case the K integral has been separated into its magnitude and angular dependences, and 
fj,'^ = cos6[, etc. Now the evaluation and simplification of the angular integral is a straightforward but lengthy 
process. The various factors of /j, may be absorbed into the Legendre polynomials using the recursion relation 



(2/ + l)xPi{x) = {l+ l)Pi+i{x) + IPi-iix). 



(7.17) 



The angular integral may be performed using the same procedure as in the scalar case: replace the Legendre poly- 
nomials with spherical harmonics using the addition theorem, and use the orthogonality of the spherical harmonics 
to eliminate the integrals. Then converting the remaining spherical harmonics back to Legendre polynomials gives 
many terms proportional to -P;(qi • (12) with various indices /. The original factors of qi • q2 in Eq. (7.16) can now be 
absorbed into the Legendre polynomials using the above recursion relation. At this point, the terms can be collected 
together, and after much algebra the formula finally simplifies to [40] 



V (/ + 2)! 
Stt (/-2)! 



dK 



A/,_2(A') 



(A-) 



(2/-l)(2/+l) (2/-l)(2/ + 3) (2/+l)(2/ + 3) 



(7.18) 



The total temperature correlation function is given by the sum of the scalar and tensor pieces as long as no correlation 
exists between the two perturbations, as will be the case for any infiationary scenario. 



B. Polarization Correlation Functions 



Direct evaluation of correlation functions involving polarization involves further complications. For correlations 
involving the polarization vector P or its magnitude, substantial simplification is not possible: since the polarization 
vector is not linear in the Stokes parameters, the x integral (replacing the ensemble average) cannot be immediately 
performed as in the temperature correlation functions above, and no progress can be made in simplifying the general 
expressions for the correlation functions. The evaluation of Eqs. (7.8) must be performed numerically through, for 
example, a Monte Carlo average over random pairs of directions separated by a fixed angle at a given point in space. 

However, the correlation functions of the Stokes parameters themselves, Eqs. (7.9), can be simplified if a small-angle 
approximation is invoked. The additional approximation is needed because, in contrast with the temperature case, 
the Q and U brightnesses have factors involving the rotation angle Consider the {QQ) correlation function, first for 
scalar perturbations. Choose the axis of the spherical coordinate system to be one of the two observation directions. 
As above, the ensemble average can be replaced by a space integration, which then eliminates one of the K integrals, 
resulting in the expression 



Q(qi)Q(q2) 



To 



To 



V 



1287r3 



5K^(2/i + l)(2/2 + l)T,,(cos(?'i)T,,(cos(?;,)A^*;^(K)A^;^(K)cos2eJ cos2e^, (7.19) 



h h 



with qi = z, which implies 6'i = 6k and (fi'i = —4>k- Now if qi and q2 point nearly in the same direction, then 
cos 2^J cos 2^2 ~ (ios24>K, using Eq. (7.3). Then the cosine factors in the above equation can be rewritten as 
(1 -|- cos4(/)if)/2. The first term then looks just like the temperature case considered earlier. The second term can be 
somewhat simplified by using the spherical harmonic summation formula to rewrite T/^ and then explicitly performing 
the integral over 4>k- The final result is 



V 

167r 



A'2rfA'|A^,(A')| 



+ -cos4,/)2Z(2/i + l)(2/2 + l) 



(/2-4)! 
ih + 4) 



A2rfAA^^^(A)A^,^ 



(A) 



(7.20) 
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A 



where the sum over I2 begins at I2 = 4. Note the final result depends explicitly on the azimuthal angle (f) since the 
value of Q changes if the axes are rotated. 

For tensor perturbations, the drill is now familiar. The assumption that Aqj(K) = Aqj(K) with the above small- 
angle approximation gives the expression 

V 



/Q(qi)Q(q2) 



\ -'0 



To To 



1287r3 



h h 



^(2/1 + l)(2/2 + 1) / K'dKA+liK)A+,^iK) 



X / dn 



K 



Pi^ (cos 6''i)T,, (cos e'2)[{l + cos^ ^'i)(l + cos^ O'^) cos^ 2(j>K + 4 cos e[ cos 0'^ sin^ 2(j>K 



(7.21) 



Then the factors of cos 9 can be absorbed into the Legendre polynomials using the recursion relation. The final 
expression for tensor perturbations is 



A'2rfA'(|5f|V4|5/| 
/o ^ 

4 fV('2+4)! '^'=jo 

'1'2 



(7.22) 



5/ (A) = (/ + 1)A+ ,(A) + /A+ ,(A) 



5f(A') 



6/^ + 9/' - /-2 /(/-I) ~ 

(2/ + 3)(2/+l) ' 2/-1 '-2 



A?t2(A'). 



The expressions for {QT) can be obtained in the same way; they are given in Ref. [43]. For gaussian initial 
perturbations, the polarization map can be reconstructed from the temperature map and the correlation functions 
{QQ) and {QT) , as shown in Ref. [20]. In principle, all of the correlation functions involving the polarization vector are 
obtainable from these as well, although there appear to be no simple formulas connecting the two sets of correlation 
functions. 



VIII. DISCUSSION 



The detection of microwave background polarization is very difficult. The most optimistic scenarios predict polar- 
ization no larger than 10% of the temperature fiuctuations, or a few parts in a million of the temperature. Needless 
to say, this sensitivity is hard to attain. The same backgrounds which affect the temperature measurements will also 
affect polarization measurements. One advantage of measuring polarization which has been realized for a long time is 
that the experiment can chop between two orthogonal polarizations on the same patch of sky, which involves rotating 
a polarizer, instead of mechanically repointing the telescope. In practice, temperature measurements can only chop at 
a few Hz, while polarization measurements can chop at hundreds of Hz. Any atmospheric noise is thus supressed much 
more effectively. A mitigating effect is that since the orientation of the horn is important, side lobes from the ground 
and diffraction effects can add noise differently to the two polarization channels. As with temperature measurements, 
the ultimate experimental limit is astrophysical foreground sources, particularly our own galaxy. 

On the theoretical side, where does polarization fit into the systematic investigation of microwave background 
anisotropies? In cosmological scenarios which invoke adiabatic initial perturbations, as in Cold or Mixed Dark 
Matter, the temperature correlation function, Eq. (7.7), contains features from which cosmological information may 
be extracted, provided reionization did not occur too early. If the universe evolved according to this broad class 
of models, measurement of the anisotropy spectrum will provide a detailed picture of the early linear growth of 
perturbations. In this case, the expected polarization can be calculated and its detection at the predicted level will 
provide important confirming evidence for the theory. Conversely, the angular scale and amplitude for the peak of 
the polarization correlation function, likely the easiest polarization quantity to measure, can serve as an additional 
piece of information when constraining theories. Also, precision calculations of the temperature anisotropies in these 



22 



models should incorporate the polarization contribution to the temperature source term in Eqs. (6.7) and (6.9); the 
polarization may affect the power spectrum at the few percent level [42]. 

If the universe did not start out with adiabatic fluctuations, as in topological defect models, or underwent early 
reionization, polarization fluctuations may take on added importance in understanding the nature of the linear pertur- 
bations, because less information will be encoded in the temperature fluctuation spectrum. Because of computational 
difficulties in the cosmic string scenario, currently the most viable defect model, no detailed temperature fluctuation 
spectra including all relevant physical processes has yet been produced, and polarization has not been investigated. 
Owing to the more complex nature of the underlying physics, it is reasonable to expect that cosmological information 
will be more difficult to recover from the cosmic string temperature fluctuation spectrum than in the adiabatic case, 
and then polarization may provide useful additional information. Likewise, if the universe underwent early reioniza- 
tion prior to z = 20, much of the information in the adiabatic microwave background spectrum will be washed out. 
In this case, the level of polarization will be signiflcantly larger [9], making it easier to detect, and may partially 
compensate for the temperature spectrum's loss of information. Additionally, in a reionization scenario, substantial 
secondary anisotropies in both temperature and polarization may be produced via the Ostriker-Vishniac effect [24]. 
Second-order effects on polarization have not yet been systematically investigated. 

Polarization limits or detections may also give useful information about intergalactic magnetic flelds at early times. 
Speciflcally, if a net circular polarization is ever detected and convincingly separated from foreground contamination, 
this would be a strong suggestion for a magnetic held polarizing the free electrons. Even upper limits on polarization 
provide one of the only methods of constraining the primordial magnetic held. 

Finally, since temperature and polarization couple differently to tensor and scalar metric perturbations, it may be 
possible to separate the two contributions with a combination of temperature and polarization measurements. Initial 
steps in this direction have been taken in Ref. [43], which considers differences in {QT) between tensor and scalar 
perturbations. The overall size of the correlations are of course small, but some differences do appear in the cross- 
correlation function. A different correlation function or some other type of statistic optimized to look for correlations 
peculiar to tensor perturbations, may give more promising results. 
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